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Abstract 



I I Using the renormalization group approach proposed by Mihis for the itinerant 

j^ ' electron systems we calculated the specific heat coefficient 7(T) for the mag- 

C^ ■ netic fluctuations with susceptibility x ~ |(5 + u;|° + f{q) near the Lifshitz 

point. The constant value obtained for a = 4/5 and the logarithmic temper- 

C I ature dependence, specific for the non-Fermi behavior, have been obtained in 

Q 



agreement with the experimental data. 



Typeset using REVT^ 



The occurrence of the non- Fermi behavior in the systems of fermions coupled to a critical 
fluctuations mode has been suggested in connection with the neutron experimentsEJB and 
studied in the framework of many body theoryScI in the case of two-dimensional (2D) and 
three dimensional (3D) models. 

Recent experiments on the heavy fermions systems showed also a non-Fermi behavior of 
these materials at low temperatures and it was associated with the proximity of quantum 
critical point (QCP). The most studied exampleu is CeCue-xAuj; where at the QCP x = 0.1 
the resistivity increases linearly with temperature T over a wide range of T and the specific 
heat C(T) is proportional to T In Tq/T. This behavior has been explainedH by the coupling 
of 3D fermionic excitations to the 2D critical ferromagnetic fluctuations near the QCP. 

The inelastic neutron scattering measurements performed on this materialsQ'Q showed the 
following new points in the behavior of this material 

1. The inelastic neutron scattering data can be fitted with a susceptibility of the form 
X"^ = C^^[f{q) + {aT — iu;)"] where a = 4/5 and not 1 as is predicted by the mean 
field approximation 

2. The quadratic stiffness vanishes, fact which shows that we are dealing with a quantum 
Lifshitz point (QLP) 

3. The peaks for x = 0.2 and x = 0.3 can be considered as 2D precursor of 3D order 

4. The scaling analysis showedEI that 'j{T) = C{T)/T has the form 

7(T) ~ T^D-l/2)a/2^1 ^-^^ 

which for D = 3 and a = 4/5 gives a temperature independent value. 

This analysis has been performed taking as the most important contribution to x the form 
containing uj" and the q-dependence of the form f{q) = Dq?^ + Cg| where a = 2/z, z being 
the critical exponent from the dynamical critical phenomenal. 

In this paper we will show that using the HertzEl renormalization group method (RNG) 
extended for T 7^ by Millialil we can obtain the In Tq/T term as a quantum correction 
to the classical results expressed by Eq. (|l]). We start from an interacting Fermi system 



and by introducing the Bose field $((?) (associated with the magnetic fluctuations) via the 
Hubbard-Stratanovich transformations one can integrate out the fermions and expand the 
effective action up to forth order in $ as 

JO /Vr JO [^TT) JO ZTT L " J 

Jo 27r io (27r)2 Jo 27i ^ ' ^ ' 

where m > is the couphng constant. The scahng variables of the models are 5" {6 is the 
deviation of the control parameter x from its critical value Xc) the stiffness D, the temperature 
and the coupling constant u. Using the transforms 

uj' = b'^'^uj q[ = hq\\ q'^ = q^Vb 

5' = h'^l'^d D' = bD (3) 

and interacting out on the shell A > g > A/6 (6 > 1) near the cut-off A we obtain the RNG 
equations. In order to calculate the speciflc heat we will use the scaling procedure to obtain 
an equation for the free energy, deflned for the free fluctuations as 



1 dz /"i d^qii /"i dq± z A sin ^ 

, — / -, — TTT / coth — - arctan — ^ . 

27r7o (27r)2 7o 2tt 2T Acos0 + gj + gj 



F = I 7^ / TTTif^- / TTT '^'-'^^ TvT^ arctan -;; — ^ — 2 r (^) 



where 6 = a tan ^ {z/5) and A " = (5^ + z'^Y''^. Following the same procedure as in Ref. 11 
we obtain the equations 

dT{b) 2 



dlnb a 



Tib) (5) 



du{b) /3 2 ^ 2 



dlnb V 2 « 
cir(6) 



_hi(6)-M^(n + <5)/2 (6) 



dlnb 

dD{b) 
dlnb 



25°(6)+2u(6)(n + 2)/i (7) 



D{b) 



''^'^-^'^l)Fib)^f. (9) 
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where n is the number of the field components, /i = fi\T{h),5°'{h),D{h)], /2 = 
/2[T(6), (5"(6), D(6)] and /a = f^\T{h)] are comphcated functions but will be approximated 
as presenting a weak dependence of 5" (6) and D{h) for 6°'{b), D{b) <^ 1. The renormalization 
procedure is stoped at 



and from Eqs. (H)-® we get 



5" (6) = 1 



T(6) = Tfe^A 



(10) 



;ii) 



u 



(b) = M63/2-2/C 



(12) 



(5"(6) = e 



2x 



(5" + 2u{n + 2) r rfxe-^(^/2+2/a) j^ (^Te^^/") 



These equations will be analyzed in two regimes. First regime defined by 



T(6) < 1 



will be called "quantum regime" and the second defined by 



T{b) > 1 



(13) 



(14) 



(15) 



will be called "classical regime". Eqs. (l^)-(lg) have been solved following Ref. |Tl| in the 
quantum and classical regimes and the renormalization coupling constant has been obtained 
as 

T^3o/4 



U 



(6) = mT^-3"/4 



l3 



(16) 



(5" + (A-4fi)(n + 2)MTi+'^/4j /2 

where 6°" = 5°'{b) and b is defined by T(b) = 1 and is T~"'^. The coupling constant have to 
satisfy the condition u{b) <^ 1 and this condition is not satisfied if 



r + (A - 45) (n + 2)mT^+"/^ = 



(17) 



If we define the coherence length by .^ ^ ~ 5 we get ^ ^ ~ j^a+i/i^ 



In order to calculate the specific heat and 7 = C{T)/T we will use Eq. (|) for the free 
energy F. Neglecting in the lower approximation the second term we obtain 

F{T) = F(6)r2/"-^/2 ^18) 



The exact solution of Eq. (|^) has the form 

/O 

and in order to perform the integral in this expression we take the variable x in the domains 



F{h) = 62/"+^/2 1'°' t/a;e-(2/"+^/2)^73 (Te^^/") (19) 



a , 1 , , 

< a; < - In - 20 

2 T ^ ' 

a 1 

-ln-<x<ln6* (21) 

2 T ^ ' 

where h* is defined as in Ref. 0, by 5{h*) = 1. In the first domainEJ fz{T) = C^T"^ and in 
the second f^{T) k, DT. Using these approximations we obtain from Eq. ([I9| ) 



^-2/a-5/2p(^jj-^ ^ ^rpl+5a/4 



1 pTb*^''" 



(22) 



where Ti = Texp [2a;/a]. If we take a = 4/5 from Eq. (^) we calculate 



F(T) = -CsT^ In - + -DT^ - -DTh*'^''^ (23) 

5 T 5 5 



and 



7(T)=7co + 7ln^ + o(^) (24) 

a result which shows that the constant value from Ref. |^ is reobtained but the RNG equations 
give also the second term specific to the non-Fermi behavior. 

Recently, Ramazashvilit3 used the same method studying QLP for such a model with 
a = 1. Our results are consistent with the results obtained in Ref. |1^ but the specific heat 
coefficient obtained has a T^^^ dependence which is given by the simple gaussian approxi- 
mation and the value a = 1. 
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